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ABSTRACT 


Tidal heating of exomoons may play a key role in their habitability, since the 
elevated temperature can melt the ice on the body even without significant solar 
radiation. The possibility of life is intensely studied on Solar System moons such 
as Europa or Enceladus, where the surface ice layer covers tidally heated water 
ocean. Tidal forces may be even stronger in extrasolar systems, depending on 
the properties of the moon and its orbit. For studying the tidally heated surface 
temperature of exomoons, we used a viscoelastic model for the first time. This 
model is more realistic than the widely used, so-called fixed Q models, because 
it takes into account the temperature dependency of the tidal heat flux, and 
the melting of the inner material. With the use of this model we introduced 
the circumplanetary Tidal Temperate Zone (TTZ), that strongly depends on the 
orbital period of the moon, and less on its radius. We compared the results 
with the hxed Q model and investigated the statistical volume of the TTZ using 
both models. We have found that the viscoelastic model predicts 2.8 times more 
exomoons in the TTZ with orbital periods between 0.1 and 3.5 days than the 
fixed Q model for plausible distributions of physical and orbital parameters. The 
viscoelastic model gives more promising results in terms of habitability, because 
the inner melting of the body moderates the surface temperature, acting like a 
thermostat. 


1. Introduction 


No exom oons have been disco vered yet, but these measurements are expected in the 


next decade. 


Bennett et af 


(120141) present a candidate, which has been detected via the 


MOA-2011-BLG-262 microlensing event. The best-£t solution for the data implies the 
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presence of an exoplanet hosting a snb-Earth mass moon. This measnrement however needs 
conhrmation, since an alternate solntion is also presented. Nevertheless, this measnrement 
indicates that the era of exomoon detections is abont to begin. 


The most favorable method for exomoon discoveries is photometry. An exoplanetary 


transit may enlighten the presence of a moon in the light curve. Details o f t 


lis method are 


thoroughly discusser 


2010 


Kipping et ah 


in the literature 


2012 


Simon et ah 


Simon et al. 


201211 . 


20071: 


Kipping 


2009a 


bl; 


Simon et al, 


I n addition, habitability of exomoons is un d er examination a s well fsee e.g. 


2010 


Heller fc Barnes 


2013i: 


Heller et al. 


201 


Hinkel fc Kanel (1201311 investigated the 


Kaltenegger 


influence of eccentric planetary orbits on moons, and concluded that a moon with sufficient 
atmospheric heat redistribution may sustain suitable temperature for life on its surface 
even if it orbits a planet that moves temporarily outside of the habitable zone (HZ) at each 
orbital period. 


Solar System analogs may serve as useful examples for different exomoon types. The 
satellites in the Solar System are diverse and life on them is a puzzling question. The 


icy surface of Europa and Ence 


a sui t able environmen t for life fjCarr et al. 


adus probably covers water ocean, which may provide 


1998 


29971: 


less et al 


Kargel et al. 


2999 


Collins fc Goodman 


291411 . Tidal and radiogenic heat keeps the interior of the body warm. 


hence maintain the water in liquid state. In fact, these internal heat sources drive to the 


erupt ion of plumes on 


well flPorco et al. 


2996 


nceladus, and s imilar phenomenon was discovered on Europa as 


Roth et al. 


291411 . 


The idea of a circumplanetary, tidally-heated habita 


inves tigated by several authors (e.g. 


Reynolds et al. 


1987 


e zone has e merged and was 


Scharl 


2996; 


Heller fc Barnes 


291311 . For the hrst time, we apply a viscoelastic model for studying tidal heat in exomoons. 
This work aims to give a detailed study of the circumplanetary Tidal Temperate Zone, and 
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discusses the differences with other models. 


Viscoelastic model 


2.1. Advantages 


Tidal heat rate of a moon is usua^ 


Reynolds et al. 

1987; 

Mever & Wisdom 

2007) 


ly ca lculated by the following expression (e.g. 




21 ko GM^R^ne 


p 


2 Q 


( 1 ) 


where G is the gravitational constant, Mp is the mass of the planet, R^, n, e and a are the 
radius, mean motion, eccentricity and semi-major axis of the moon, respectively. Q is the 
tidal dissipation factor and k 2 is the second order Love number: 


ko = 


3/2 


1 + 


19m 

2pgRn 


( 2 ) 


where fi is the rigidity, p is the density and g is the surface gravity of the satellite. This 
calculation method is called the fixed Q model, because Q, p and k 2 are considered to be 
constants. 


The hxed Q model is broadly used in ti c 


the tidal heat of the body (IRoss fc Schubert 


al ca l culations, but highly un derestimates 


1988 


Mever fc Wisdom 


2007). Moreover, 


both Q and p are very difficult to determine, and vary on a large scale for different 
bodies: from a few to hundreds for rocky p lanets, and tens or hundreds of thousands 


for giants (see e.g. 


Goldreich fc Soter 


19661) . In addition, these parameters are not 


constants, since they stron . 


2003 


Henning et al. 


2009 


glv depend on t 

le temperature ( 

Fischer & Soohn 

1990; 

Moore 

Shoii & KnritJ 

2014 

). As a consequence, tidal heat flux has 
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a temperature dependency, as well; it reaches a maximum at a critical temperature (Tc) 
as can be seen in Fig. [TJ Between the solidus and the liquidus temperature (Tg and Ti, 
respectively) the material partially melts. Above the breakdown temperature (Tb) the 
mixture behaves as a suspension of particles. The dashed curve represents the convective 
heat loss of the body. Circles indicate equilibria, for example, the solid circle between Tc 
and Tb is a stable equilibrium point. If the temperature increases, convective cooling will 
be stronger than the heat flux, resulting in a cooler temperature. In case of decreasing 
temperature, the tidal heat flux will be the stronger, hence the temperature increases, 
returning the system to the stable point. The stable equilibrium between the tidal heat and 


convection is not necesserely located between Tc and Th, in fact, t 


two curves do not have intersection at all (see 


Henning et ah 


lere are cases, when the 


20091. Fig. 6). In these cases 


tidal heat is not strong enough to induce convection inside the body. 


In contrast to the hxed Q model, viscoelastic models take into account the temperature 
dependency of the body, hence are more realistic. 


2.2. Description 


In viscoelastic models k 2 /Q is replaced by the imaginary part of t he complex 
number Im(/c 2 ), which describes structure and rheology in the satellite fjSegatz et af 


jOve 


1988h 


21 R® 


( 3 ) 


Note that in this expression the mass of the planet and the semi-major axis of the moon 
are eliminated by the mean motion (n = ^jGMp/a^). 


Henning et ah 


020091 ) gives the value of Im(fc 2 ) for four different models (see Table 1. 


in their paper). In this work we use the Maxwell model; 
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- /m(/c2) 


57rj0J 


ApgR^ 


1+1 + 


19 m 

‘IpgR-a 


(4) 


where rj is the viscosity, oj is the orbital frequency and p is the shear modulus of the 
satellite. 


The viscosity and the shear modulus of the body strongly depend on the temperature. 
Below the Tg the shear modulus is constant: p = 50 GPa and the viscosity follows an 
exponential function: 


r] = rio exp 


E 


(5) 


where r^o = 1.6 • lO^Pas (reference viscosity), E is the activation energy, is t he universal 


gas constant and T is the temperature of the material flFischer fc Spohn 


199nh . 


Between Tg and Tb the body starts to melt. The shear modulus changes by 


where = 8.2 • 10‘^K and /i 2 
expressed by 


^ = 10(T+")Pa, 


( 6 ) 


—40.6 flFischer fc Spohn 


19901). The viscosity can be 


?7 = ? 7 o exp ( ^ j exp {-Bcf)) , 


(7) 


where 0 is the melt fraction which increases linearly with the tempe rature 


7] (0 < 0 < 1) and B is the melt fraction coefficient (10 < B < 40) ( Moore 


Detween Tg and 


20031). 


At Tb the grains disaggregate, leading to a sudden drop in both the shear modulus 
and the viscosity. Above this temperature the shear modulus is set to a constant value: 
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/i = 10 ^Pa. The viscosity follows th e Roscoe-Ein stein relationship so long as it reaches the 
liqnidns temperatnre (where 0 = 1) flMoordl2003l) : 


/40000K\ , ,_5/2 

rj = 10 ^exp ( ——— j (1.350 — 0.35) Pas . 


( 8 ) 


V T 

Above Ti the shear modulus stays at 10“^Pa, and the viscosity is described by (iMoore 


200311 


,= 10-Vp(^2^)pa = . 

In our calculations rocky bodies a re considered a s sate llites, and for this reason 


( 9 ) 


we follow the melting temperatures of 


Henning et ah 


(2009), namely: Tg = 1600 K, 


Ti = 2000 K. We assume that disaggregation occurs at 50% melt fraction, hence the 
breakdown temperature will be Tf, = 1800 K. 


2.3. Internal structure and convection 


The structure of the moon in the model is the following: the body consists of an inner, 
homogenous part, which is convective, and an outer, conductive layer. If the tidal forces are 
weak, the induced temperature will be low, resulting in a smaller convective region and a 
deeper conductive layer. But in case of strong tidal forces, the temperature will be higher, 
hence the convective zone will be larger with a thinner conductive layer. 


For calcul 

ating t 

Henning et ah 

(2009) 


Qbl 


= k 


- mantle 




surf 


therm" 


S{T) 


( 10 ) 





















where /ctherm is the thermal conductivity (~ 2W/mK), Tmantie and Tgurf are the temperature 
in the mantle and on the surface, respectively, and S{T) is the thickness of the conductive 
layer. We use S{T) = 30 km as a hrst approximation, and then for the iteration 


6{T) 


_d_ f Ra\ 
202 \RacJ 


( 11 ) 


is used, where d is the mantle thickness (~ 3000 km), 02 is the flow geometry constant 
(~ 1), Rqc is the critical Rayleigh number (~ 1100) and Ra is the Rayleigh number which 
can be expressed by 


Ra = 


agpd'^ gsL 


( 12 ) 


f]iT') K kthevm 

10“'^) and K is the thermal diffusivity: n = /cthprm/fpC'p) 


Here a is the thermal expansivity 
with C n = 1260 J/(kg K). For detailed description see the clear explanation of 

mm . 


Henning et ah 


Because of the viscosity of the material the thickness of the boundary layer and 
convection in the underlying zone changes strongly with temperature. The weaker 
temperature dependencies of density and thermal expansivity are neglected in the 
calculations. The iteration of the convective heat flux lasts until the difference of the last 
two values is higher than 10“^°W/m^. 


Calculations of tidal heat flux and convection are made for a hxed radius, density. 


eccentricity and orbita’ 
the equilibrium state. 


period of the moon. We assume that with time, the moon reaches 


Henning et ah 


(120091) showed that planets with significant tidal 


heating reach equilibrium with convection in a few million years. However, change in the 
eccentricity can shift, or destroy stable equilibria. After finding the stable equilibrium 
temperature, the tidal heat flux is calculated, from which the surface temperature can be 
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obtained using the Stefan-Boltzmann law: 


^surf — 




1/4 


tidal 


ATiR'la 


(13) 


where a is the Stefan-Boltzmann constant. This is the hrst time of using a viscoelastic 
model for obtaining the tidal heat induced surface temperature on exomoons. 


2.4. Results 


The satellite’s surface temperature is calculated for different orbital periods and 
radii, at a hxed density and eccentricity. Stellar radiation and other heat sources are not 
considered, and have been neglected. The orbital period and the radius of the moon varies 
between 2 and 20 days, and between 250 km and 6550 km, respectively. It is common to 
consider Earth-mass moons in extrasolar systems when speaking of habitability, however, 
their existence is not proven. In the Solar System the largest moon, Ganymede has only 


0.025 E arth mass. But 
planet. ICanup fc Ward 


he ma ss of satellite systems is proportional to the mass of their host 


(120061) showed that this might be the case for extrasolar satellite 
systems as well, giving an upper limit for the mass ratio at around 10“^. This means, that 
10 jupiter-mass planets may have Earth-mass satellites. Besides accretion, large moons 
can also form from collisions, as in the case of the Earth’s Moon. Other possibility is the 
capturing of ter restrial-sized bodies through a close planetary encounter, as described by 


WilliamsI (120131) . For these reasons, we also take Earth-like moons into account. 


The results can be seen in Fig. |2l where the density of the moon is that of lo, and its 
eccentricity is set to 0.1. Different colours indicate different surface temperatures. In the 
white region there is no stable equilibrium between tidal heat and convective cooling. In 
other words, tidal heat is not strong enough to induce convection. For comparison, a few 
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Fig. 1.— Schematic fig ure of the temperature dependency of tidal heat flux and convective 


heat loss flMoore 


2003|). 
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2.— Tidal heat induced surface temperature for moons with similar density to lo at 0.1 


orbital eccentricity. 



















11 


Solar System moons are plotted that have similar densities to lo’s. Yellow contonr cnrves 
denote 0 and 100 °C. The green area between these cnrves indicates that water may be 
liqnid on the snrface of the moon (atmospheric considerations were not applied). We define 
this territory as the Tidal Temperate Zone (TTZ). 

Interestingly, the location of the TTZ strongly depends on the orbital period, and less 
on the radins of the moon. Low radii are less relevant, since smaller bodies are less capable 
of maintaning signihcant atmospheres. 

The dependency on the eccentricity can be seen by comparing Fig. [2] and El In the case 
of the latter fignre the moon’s eccentricity is 0.01. For most of the orbital period-radius 
pairs there is no solntion (white area). Dne to this drastic difference, Enropa analognes 
get ont of eqnilibrinm for smaller eccentricities, and the TTZ becomes narrower and shifts 
to shorter orbital periods. Note, that radiogenic heat is not considered in the model, 
which conld pnsh back the moon into eqnilibrinm state, and wonld resnlt in higher snrface 
temperatnre. 

Similar calcnlations were made for the density of the Earth and Titan (left and right 
panel of Fig. EJ respectively). The densities do not have high inflnence on the tidally 
induced snrface temperatnre, however, the TTZ slightly shifts to lower orbital parameters 
for higher densities. (The density of Earth, lo and Titan are 5515 kg/m^, 3528 kg/m^ and 
1880 kg/m^, respectively.) 

In the left panel of Fig. Elan example Earth-like moon is plotted inside the TTZ. This 
hypothetic body has the same mean snrface temperatnre (288 K), radins (6370 km) and 
density as the Earth, hence its orbital period is 2.06 days. In the right panel a few Solar 
System satellites are plotted that have similar densities to that of Titan. 

The stellar flnx for moons with ambient temperatnres of ~100 K (which is similar to 
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Fig. 3.— Tidal heat induced surface temperature for moons with similar density to lo at 
0.01 orbital eccentricity. 
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the case of the Galilean and Satnrnian moons in the Solar System) is abont one percent of 
the tidal flnx in the TTZ. For this reason, stellar insolation may be safely ignored if the 
planet-moon system orbits the star at a far distance, or if they are free-floating. For moons 
system in which the stellar irradiation alone is snfficient to heat the snrface to levels of 
order the melting temperatnre or higher, the models presented here are wonld need to be 
replaced by more complex hybrid ones to take into acconnt both sonrces of heat and their 
very different spatial distribntions on and within the moon. 


3. Comparison to the fixed Q model 
3.1. Method 


It is clear from the results, that the viscoelastic model does not give solution in case of 
small tidal forces. In other words, the amount of heat that is produced by tidal interactions 
is insufficient to induce convective movements inside the body, and for this reason there 
is no equilibrium between them. In contrast, the hxed Q model provides solution in these 
cases, as well. However, the viscoelastic model describes the tidal heating of the body more 
realistically than the fixed Q model, due to the temperature dependency of the Q and fi 
parameters. How are the results of the two models related to each other? 


from 


or comparing th e resu 


ts of the two kinds of models, we use the expression of Eq. (7) 
Peters fc Turner! fl2013l) for the hxed Q calculation: 


9747^2 ) \ liQ ) V/915/2^ j ’ 


( 14 ) 


where Tgurf is the surface temperature of the moon induced by tidal heating, G is the 
gravitational constant, a is the Stefan-Boltzmann constant, is the radius, p is the 
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densitiy, /i is the elastic rigidity, and Q is the dissipation fnnction of the moon, e is the 
eccentricity of the moon’s orbit, and [3 is expressed with the semi-major axis (a) and the 
mass of the planet (Mp): 

/Q/bf \ 

a = Pan = P j , (15) 

where is the Roche radius of the host planet. These equations can be used for calculating 
the surface temperature of the moon heated solely by tidal forces. 

The viscoelastic model is described in details in Sections 2.2 and 2.3. The satellite’s 
mean motion can be expressed from /3 by 


2tiG p 


which makes the comparison of the two models easier. 


(16) 


3.2. Surface temperature 


For comparison of the hxed Q and the viscoelastic model, see Figs. [5l [6]and[71 which 
show the surface temperature of a moon, calculated with both the viscoelastic (red solid 
curve) and the hxed Q model (green dashed curve) as functions of the eccentricity, radius 
and orbital period of the satellite. For the density of the moon we used 5515 kg/m^, 
which is the density of the Earth, a nd for the hxed Q mo del we used Q = 280 and 


p = 12 ■ 10^°kg/(ms^) in each case (iPeters fc Turner! 120131 . Table 1). The radius, orbital 


period and eccentricity of the satellite are set to that of the Earth, lo and 0.03, respectively, 
except that one of these parameters is varied in each hgures (horizontal axes). The 
horizontal light blue, dashed lines indicate 0 and 100 °C (making the boundaries of the 
TTZ), and the solid blue lines denote the minimum and maximum temperatures (—20 and 
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60 °C ) that are probable limits of habitability on an Earth-like body flSnllivan fc Baross 


20071. Chapter 4). In salty solntions the lower limit for microbial activity is aronnd —20°C, 


and the npper limit for complex eukaryotic life is 60 °C. The latter temperature is also 
about the runaway greenhouse limit for Earth. These limits are only used for Earth-like 
bodies (p = pEarth and ~ .REarth)- Vertical lines show a few examples from the Solar 
System for different eccentricities, radii and orbital periods. 


It is noticeable that the red curve is less steep than the green curve, and larger portion 
of it is located inside the TTZ, especially in Figs. [5] and [Hi It shows that the viscoelastic 
model stabilizes the surface temperature comparing to the hxed Q model. These are just 
a few examples indicating that the viscoelastic model is less sensitive to these parameters, 
and that there are huge differences in the results of the models. In the next section the 
volume of the TTZ is investigated more thoroughly. 


3.3. Occurrence rate of ’habitable’ moons 

Habitability on extraterrestrial bodies is an exciting, but complex question. Here we 
consider solely the tidally induced surface temperature of a hypothetic moon. We were 
curious about the occurrence rate of moons with suitable surface temperature for life. For 
this reason we mapped the phase space evenly with hypothetic moons that have different 
radii (between 250 and 6550 km) and eccentricities (between 0.001 and 0.1), and their 
densities are that of the Earth. We used both the viscoelastic and the hxed Q model for 
calculating the surface temperature of these bodies, and then calculated the percentage 
of those that have suitable surface temperature, i.e. that are located inside the TTZ 
(0 < Znrf < 100 °C). The calculation was made for different orbital periods between 0.1 
to 3.5 days, and for each value there were 63100 hypothetical moons distributed in the 
radius-eccentricity phase space. The result can be seen in Fig. |8l Red solid and green 
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Fig. 5.— Tidally induced surface temperature of the satellite as fuction of its eccentricity. 



Fig. 6.— Tidally induced surface temperature of the satellite as fuction of its radius. 
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dashed curves indicate the percentage of being inside the TTZ for the viscoelastic and for 
the hxed Q model, respectively. The blue dotted curve shows the percentage of those cases 
that do not give result for the viscoelastic model. The top axis shows the /3 parameter, 
which is the ratio of the moon’s semi-major axis and the planet’s Roche-radius. It can be 
clearly seen that the red and green curves have a peak, which means that the probability 
of having suitable surface temperature has a maximum at a certain orbital period. The 
viscoelastic model predicts a much more efficient heating than the hxed Q model, i.e. a 
much larger fraction of the hypothetical moons have their surface temperature between 0 
and 100 °C. The ratio of the integral under the red to that under the green curve is 2.8, 
meaning that 2.8 times more exomoons are predicted in the TTZ with the viscoelastic 
model. For the viscoelastic model the maximum percentage appears around 1 day orbital 
period, and here the probability for the moon of being inside the TTZ is almost 80%. For 
higher orbital periods, this probability rapidly falls down, which is in contrast with the 
hxed Q model. The latter has its peak around 1.5 days, and it has less than 20% chance 
for satellites being in the TTZ. Despite of the high probabilities achieved by the viscoelastic 
model for small orbital periods, the hxed Q model give more promising results for those 
moons that have their orbital periods at 2 days or more. 

For detailed study, the 0 and 100 °C temperature contours were plotted in the 
radius-eccentricity plane for a few, specihc orbital periods, namely P = 0.5 day (top panel), 
P = 1 day (middle panel) and P = 1.5 days (bottom panel) (see Fig. [9]). Again, red and 
green colours represent the viscoelastic and the hxed Q model, respectively. Between the 
contour curves the region of the TTZ is hlled with light red and light green colours. The 
result shows that the viscoelastic model mostly favours the small moons, especially at high 
eccentricities, but also some large moons at small eccentricities over the hxed Q model. 
This suggests that the viscoelastic model is less sensitive to the parameters of the moon, 
and holds the temperature more steady than the hxed Q model. This is due to the melting 
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of the inner material o 


the moon that leads to a less elevated temperature, as discussed by 


Peters fc Turner! (120131) . On the other hand, the lower temperature implies that the total 


irradiated flux of the moon will be also lower, hence making the detection of the moon more 
difficult. 


We were also interested in ’Earth-twins’ as satellites and in the probability of their 
’habitability’. For this reason we made similar calculations, but the radius and the density 
of the hypothetical moons were set to be close to that of the Earth: = 6378km(±5%) 

and p = 5515kg/m^(±5%). The radius and density values were chosen randomly from 
these intervals. The eccentricity was altered similarly as in the previous case (uniformly 
between 0.001 and 0.1). Altogether 200000 cases were considered for each orbital period. 
The temperature limits were set to —20 and 60 °C, which are the probable limits for life 
on Earth. Fig. [TOl shows the results of this calculation. Note, that the peaks of the red 
solid (viscoelastic model) and the green dashed (fixed Q model) curves are shifted to higher 
orbital periods, comparing to Fig. [HI This is in part caused by the changed temperature 
limits, and in part by the much shorter radius range. The maximum probabilities are also 
higher, that is especially visible in the case of the fixed Q model, which reaches more than 
40 % at the curve’s peak (in the previous case it was less that 20 %). As one would expect, 
it suggests that larger moons are more probable of maintaining warm surfaces. The ratio of 
the areas under the red and the green curves is 2.3. 

In general, it can be concluded that the viscoelastic model is not just more realistic 
than the fixed Q model, but also gives more promising results for exomoons, since much 
larger fraction of the hypothetical satellites have been found in the TTZ. In those cases 
when the viscoelastic model does not give solution for the equilibrium temperature, one can 
use the fixed Q model instead, however, the values of Q and p are highly uncertain. 
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3.4. The value of Qfi 


With the product of Q and /i, one can easily calculate the tidally induced surface 
temperature of a moon without using a complex viscoelastic model. Using Eq.[Tl] is a 
fast way to obtain Tg^ri, but a good approximation is needed for the Qfi value. For such 
calculations, in the following, we give the Q/i values for hypothetical moons. Because of the 
large number of possible variations in the physical and orbital parameters of the moons, 
only a few. Solar System-like bodies were considered. Since the Qn varies several orders of 
magnitudes for different rocky bodies, a good estimation can serve almost as well as the 
exact value. The following examples can be used as a guideline for making such estimations. 
Note, that the used model can be applied to rocky bodies, but for icy satellites, such as 
Enceladus or Europa, the results may be misleading, because of the more complex structure 
and different behaviour of icy material. 


From the surface temperature of the moon, that was calculated from the viscoelastic 
model, the Q/r value was determined using Eq.[Tl]for six orbital period-eccentricity pairs. 
In Table [T] the tidally induced surface temperatures and the logarithm of the Qfi values 
can be seen. The radius and density of the satellites are thos e of the corresponding Sola r 


System bodies (see the first column), and the values are from 


Murray fc Dermott 


(I1999I). 


The eccentricities are set to 0.01 and 0.1, and the orbital periods to 1, 2 and 3 days. ’N.a.’ 
indicates that there was no solution (weak tidal forces). 


3.5. Scaling the Galilean satellite system 

Since no satellite has been discovered so far outside the Solar System, we used the 
Galilean system as a prototype for realistic calculations, lo, Europa and Ganymede are 
orbiting in a 1:2:4 mean motion resonance, that maintains their eccentricities, which play an 
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label 

Rm [km] 

p [kg/m^j 

e 

P [days] 

Tsurf [K] 

logio(Q/i) [Pa] 





1 

281 

14.0 




0.01 

2 

213 

13.0 

Earth-like 

6378 

5515 


3 

180 

12.4 


1 

378 

15.5 




0.1 

2 

291 

14.4 





3 

249 

13.8 





1 

256 

12.5 




0.01 

2 

194 

11.5 

Mars-like 

3394 

3933 


3 

163 

10.9 


1 

342 

14.0 




0.1 

2 

263 

13.0 





3 

225 

12.4 





1 

232 

11.1 




0.01 

2 

176 

10.1 

Moon-like 

1738 

3342 


3 

n.a. 

n.a. 


1 

313 

12.6 




0.1 

2 

240 

11.5 





3 

205 

10.9 





1 

235 

11.2 




0.01 

2 

177 

10.2 

lo-like 

1821 

3532 


3 

n.a. 

n.a. 


1 

315 

12.7 




0.1 

2 

242 

11.7 





3 

207 

11.1 


Table 1: Qfi values for Solar System-like, rocky moons. The radii and d ensities are those of 


the co rresponding Solar System bodies. The reference for these valnes is 


Murray fc Dermott 


1119991. Appendix A). 
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essential role in forcing continuously their tidal heating. 


Qgihara fc Idal (120 12h investigated 


satellite formation in the circumplanetary disc of giant planets using A^-body simulations 
including gravitational interactions with the circumplanetary gas disc. They have found 
that 2:1 mean motion resonances are almost inevitable in Galilean-like satellite systems, 
and based on their results they predict that mean motion resonances may be common 
in exoplanetary systems. For these reasons the Galilean satellite system seems to be a 
representative example for realistic calculations, since the moons are in resonance, and their 
scaled-up versions will probably stay in resonance, too. 


The test systems consist of a planet (Jupiter) and the four Galilean moons. 91 
cases are considered, one is the real Galilean system, and the others are the scaled-up 
versions: the masses of the planet and the moons were multiplied by the scale factor 
{scale = 1.0,1.1,1.2,... 10.0), and the semimajor axes of the moons were altered with 
constant orbital periods for each scale value. 


scale ■G{Mp + M,) ’ 

where a is the semimajor axis of the moon, Mp and Mg are the masses of the planet and the 
moon, respectively. The hxed orbital periods guarantee that the satellites approximately 
stay in resonances. This calculation resulted in constant (3 values for all scale parameters. 

Using both the hxed Q and the viscoelastic models, the warmth of tidal heat was 
investigated in each case. The tidal heat induced surface temperature can be seen in 
Fig. [m where the 91 cases are connected with solid curves for each satellite. 182 
other cases were calculated, too, which are shown with dashed and dotted curves in the 
hgure. These curves indicate that the densities (dashed curve) and the eccentricities 
(dotted curves) of the satellites are doubled compared to their original values in the 
Solar System. In the calculations ja and Q were set for all satellites to that of lo, namely 
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10^°kg/(ms^) and 36, respectively, except for Europa, wh ich has the following parameters: 
M = 4 ■ 10^ kg/(m s^) and Q = 100 fjPeters fc Turneiil2013l) . Densities of the moons are from 


Lodders &: Fedevl (119981). an c 


valnes is 


Murray fc Dermott 


the reference for the semi-major axis, eccentricity and mass 


lOOOl. Appendix A). 


For lo, in the scale = 1 (Solar System) case the fixed Q and the viscoelastic models 
give 60 K and 160 K, respectively. The observe d surface heat flux induced by tidal heat 


on lo is around 2W/m^, which is a lower limit fjSpencer et af 


20001) ■ In other words, tidal 


forces produce at least 77 K heat on the surface of lo. The fixed Q model resulted in a lower 
value than this limit, but note that Q and are very difficult to estimate. The viscoelastic 
model gave much higher temperature than the observation, but keep in mind, that the 


heat is concentrated in hotspots, and is not evenly distributed on the surface of 


temperature of the most warm volcano, Loki is higher than 300 K (jSpencer et al. 


o. The 


20001 ) 


The viscoelastic model gives solution only for lo (orange curves) and Europa (light blue 
curves), but not for all scale values, as shown in Fig. HH In those cases when the densities 
are the twice than those in the Solar System (dashed curyes), the surface temperatures 
are just slightly higher. In fact, in the viscoelastic model, higher densities result in less 
tidal heat because of the imaginary part of the second order Love number. Doubling the 
eccentricity instead of the density (dotted curves) makes the surface temperature higher in 
each case. 


4. Conclusions 

We have used, for the hrst time, a viscoelastic model for calculating the surface 
temperature of tidally heated exomoons. The viscoelastic model gives more reliable results 
than the widely used hxed Q model, because it takes into account that the tidal dissipation 
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factor (Q) and rigidity (ju) strongly depend on the temperature. Besides, these values are 
poorly known even for planets and satellites in the Solar System. Using the viscoelastic 
model for exomoons helps to get a more realistic estimation of their surface temperature, 
and to determine a circumplanetary region, where liquid water may exist on them. It may 
help future missions in selecting targets for exomoon detections. 

We have dehned the Tidal Temperate Zone, which is the region around a planet where 
the surface temperature of the satellite is between 0 °C and 100 °C. No sources of heat were 
considered other than tidal forces. Assuming, that the planet-moon system orbits the star 
at a far distance, or the stellar radiation is low due to the spectral type, tidal heat can be 
the dominant heat source affecting the satellite. We have investigated such systems, and 
found that the TTZ strongly depends on the orbital period, and less on the radius of the 
moon. For higher densities or eccentricities of the moon, the location of the TTZ is slightly 
closer to the planet. 

Comparing this model to the traditionally used hxed Q model revealed that there 
are huge differences in the results. Generally, the viscoelastic model is less sensitive to 
moon radius than the uniform Q model, keeping the surface temperature of the body more 
steady. The reason is that higher tidal forces induce higher melt fraction which results in a 
lower temperature than the hxed Q model. The viscoelastic model demonstrates the way in 
which partially melting of a moon can act as a thermostat and tend to hx its temperature 
somewhere near its melting point over a wide range of physical and orbital parameters. 

As a consequence, the statistic volume of the TTZ is much larger in the viscoelastic case, 
which is favourable for life. But this lower temperature also means that the detectability of 
such moons is lower in the infrared. In addition, for low tidal forces there is no equilibrium 
with the convective cooling; hence, only the hxed Q model provides solution. In these cases 
the challenge is to determine the values of Q and /i. 
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For a few characteristic cases the product of the tidal dissipation factor and rigidity 
was calculated from the viscoelastic model, in order to help in quick estimations of tidally 
heated exomoon surface temperatures. Since the viscoelastic model is more realistic 
because of the inner melting and the temperature dependency of the parameters, but the 
hxed Q model is easier to use, these Qfi values (along with the surface temperature) are 
provided in Tabled] By inserting Q/i into Eq.dH one can get the estimation of the tidally 
induced surface temperature of a moon. Connection between the quality factor [Q) and 
the vis coelastic parameters (viscosity and shear modulus) was given for the Maxwell model. 


too, by 


Remus et af 


( 1201211 . 


Earth-like bodies were also investigated as satellites, and in these cases the —20 and 
60 °C temperatures were used as limits of habitability. The results are similar, but the 
volume of this habitable zone is larger than that of the TTZ for wide range of satellite radii. 
This habitable zone includes atmospheric considerations of the moon, but stellar radiation 
was neglected in the calculations. In case of signihcant radiation from other sources, the 
surface temperature of the moon will be higher. Additional heat sources (such as stellar 
radiation, radiogenic processes, reflected stellar and emitted thermal radiation from the 
pl anet), and the effects of eclipses, or the obliquity of the satellite are thoroughly discussed 


by 


Heller fc BarnesI ( 201 


3 )- 


For simulating realistic systems, the Galilean moons were used as a prototype. Their 
surface temperature was calculated with both models for different, scaled up masses. The 
mean motion resonance between the satellites helps to maintain their eccentricity, and 
consequently to maintain the tidal forces. By raising their masses, the temperatures of lo 
and Europa elevates less drastically in the viscoelastic model, than in the hxed Q model 
(see Fig. m- At scale = 5 (masses are hve times as in the Solar System case) the surface 
temperature of Europa is ~ 150 K calculated from the viscoelastic model. Assuming that 
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its density does not change, its radius will be approximately 0.25 Earth radii. In case of 
an additional 100-120 K heat (e. g. from stellar radiation), the ice would melt, and this 
super-Europa would become an ’ocean moon’, covered entirely by global water ocean. The 
used viscoelastic model might not be adequate, and can be oversimplified for such bodies 
that consist of rocky and icy layers, as well. Salty ice mixtures may also modify the results. 
The applied model ignores the structure, pressure and other effects, and applies melting for 
the whole body. However, it provides a global picture of the tidally heated moon. Even 


with a more detailed viscoelastic moc 
(rocky core, ocean and ice shell). 


, tha t describes Enceladus as a three layered body 


Barn (120081) have found that tidal heat is ~ 10 tim es lower 


than 


hat was observed by the Cassini Composite Infrared Spectrometer. Similarly, 


Moore 


(120031) concluded that observed heat flux on lo is about an order of magnitude higher than 


that can be explained with a multilayered, viscoelastic model. These results suggest that 
tidal heat can be much more relevant than what is predicted by models. 


We thank Amy Barr, Rene Heller, Edwin Kite and Omer Bromberg for useful 
conversations. We are also grateful to Laszlo L. Kiss and Darren Williams for useful 
comments that greatly improved the manuscript. This research has been supported in part 
by the World Premier International Research Center Initiative, MEXT, Japan. VD has 
been supported by the Hungarian OTKA Grants K83790, K104607 and the LenduIet-2009 
Young Researchers Program of the Hungarian Academy of Sciences, and the ESA PECS 
Contract No. 4000110889/14/NL/NDe. 
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Fig. 7.— Tidally induced surface temperature of the satellite as fuction of its orbital period. 
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Fig. 8.— Percentage of cases that give surface temperatures between 0 and 100 °C. Solid 
red and dashed green curves represent the results of the viscoelastic and the hxed Q model, 
respectively. The density of the moon is that of the Earth in each case. The dotted blue 
curve shows the ratio in percentage of those cases where the viscoelastic model did not have 
solution for the equilibrium temperature. 
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Fig. 9.— Temperature contours for the two kinds of models. Red colour: viscoelastic model, 
green; fixed Q model. Top panel: orbital period P = 0.5 day, middle: P = 1 day, bottom: 
P = 1.5 days. 
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Fig. 10.— Percentage of cases that give surface temperatures between —20 and 60 °C. Solid 
red and dashed green curves represent the results of the viscoelastic and the fixed Q model, 
respectively. The density of the moon is that of the Earth in each case. The dotted blue 
curve shows the ratio in percentage of those cases where the viscoelastic model did not have 
solution for the equilibrium temperature. 
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Fig. 11.— Surface temperature of the Galilean moons as function of the scale parameter, 
calculated by both the hxed Q and the viscoelastic model. Dashed and dotted curves indicate 
that the density or the eccentricity of the moon is doubled compared to the Solar System 
case, respectively. 
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